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ABSTRACT 

We have mapped the central 10'xlO' of the dense core TMC-1C at 450 and 
850 /im using SCUBA on the James Clerk Maxwell Telescope. The unusally 
high quality of the 450 /zm map allows us to make a detailed analysis of the 
temperature and column density profiles of the core. We find that the dust 
temperature at the center of TMC-1C is ~ 7 K, rising to ~ 11 K at the edges. 
We discuss the possibility and effects of a variable emissivity spectral index on 
the derived mass profile. The low dust temperature of TMC-1C results in a high 
derived mass for the core, significantly larger than the virial mass estimated from 
the linewidth of the N2H + (l-0) transition. This result is valid within a wide range 
of dust properties and ellipticities of the core. The N2H + (l-0) spectra, taken with 
the IRAM 30m telescope, show signs of self-absorption, which provide evidence 
of sub-sonic infall motions. The derived density profile and infall velocity is 
compared to the predictions of several popular star formation models, and the 
Bonnor-Ebert model is the best fit analytic model. 

Subject headings: stars: formation — dust, extinction — submillimeter 



1. Introduction 

Prestellar cores are self-gravitating condensations in molecular clouds that are much 
denser and colder than their surrounding medium. These objects (such as L1544 and L1544F) 
are believed to be on the verge of collapse, and infall motions are not uncommon (Williams et 
al. 1999; Crapsi et al. 2004). Because prestellar cores are at a critical point on the way to be- 
coming a star, their properties, such as density and temperature distribution, as well as their 
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velocity field, are of particular interest. The density and temperature profiles for prestellar 
cores are difficult to obtain from molecular line data, because gases can deplete onto the 
dust grains at the high (n > 10 4 cm -3 ) densities and low temperatures (T < 10 K) present 
in the interiors of cores. A more accurate estimate can be obtained from measurements of 
the dust emission, which peaks at submillimeter wavelengths and is often optically thin. Ob- 
servations at multiple submillimeter wavelengths can be used to determine such properties 
as the temperature, density, emissivity, and spectral index of the dust (Hildebrand 1983). 
Kinematic information, such as infall, outflow, and rotation, can be determined by molec- 
ular line observations of carefully chosen high and low density tracers (Caselli et al. 2002b; 
Belloche et al. 2002; Di Francesco et al. 2001). 

TMC-1C is a starless core in the Taurus molecular cloud complex, at a distance of 140 
pc (Kenyon, Dobrzycka, & Hartmann 1994). Taurus is known to be a site of low mass, 
somewhat isolated star formation. TMC-1C has previously been observed to be a "coherent 
core," meaning that its velocity dispersion is constant, at slightly more than the sound speed, 
over a radius of 0.1 pc (Barranco & Goodman 1998; Goodman et al. 1998). In addition, it 
shows a velocity gradient consistent with solid body rotation, at 0.3 km s _1 pc -1 (Goodman 
et al. 1993). Because TMC-1C is fairly nearby and has interesting kinematics, it is a good 
candidate starless core for further study. 

In this paper we use continuum data taken with SCUBA at 450 and 850 /zm to determine 
the temperature and mass distribution of the dust in TMC-1C, and spectral line maps from 
the IRAM 30m dish to provide velocity information. In Section 5.4 we note how changes in 
the assumed dust emissivity alter the conclusions of our analysis. In Section 5.6 we compare 
the mass derived from our multiwavelength observations to the mass we would have derived 
if only 850 jum data had been available to us. In Section 5.7, various collapse models for 
star formation are considered, and their predictions for infall speed and density profile are 
compared to observations. Finally, we conclude that for a wide range of assumptions, TMC- 
1C is a collapsing coherent core that is not especially well fit by any of the star formation 
models that we available. 



2. Observations 
2.1. Continuum 

We observed a 10'xlO' region of TMC-1C in Taurus with SCUBA (Holland et al. 1999) 
on the JCMT in exceptionally stable grade 1 weather. We used the standard scan-mapping 
mode, creating 850 /iin and 450 fim maps simultaneously (Pierce-Price et al. 2000; Bianchi 
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ct al. 2000). Three chop throw lengths of 30", 44", and 68" were used in both the right 
ascension and decination directions. The JCMT has beam widths of 7.5" at 450 /im and 
14" at 850 /im which subtend diameters of 0.005 and 0.01 pc, respectively, at the distance 
of Taurus. Pointing during the observations was typically good to 3" or better. 

The data were reduced using the SCUBA User Reduction Facility (SURF) (Jenness 
& Lightfoot 1998). The data were flatfielded, extinction corrected using skydips, despiked, 
baseline corrected, and had sky noise removed. Maps were then made using Emerson Fourier 
deconvolution (Emerson 1995). The resultant maps were calibrated with the source CRL618, 
with uncertainties in the calibrator's flux densities of ~3.7% at 850 /im and ~12.5% at 450 
/im. The noise level of the 850 /im map is 0.009 Jy/beam, and 0.053 Jy/beam in the 450 /im 
map. Not that the scan-mapping mode is insensitive to structures on scales greater than a 
few times the largest chop throw. The flux maps can be seen in Figure 1. 

2.2. Spectral Line 

The 3 mm N2H + (l-0) observations shown in this paper were taken with the IRAM 30m 
telescope, and will be discussed in more detail in an upcoming paper (Schnee et al. 2005). A 
beam sampled map of the inner ~2' of TMC-1C was obtained in frequency switching mode 
with a spectral resolution of 0.06 km s" 1 . The data were reduced using the CLASS package 
developed jointly between Observatoire de Grenoble and IRAM. Second order polynomial 
baselines were subtracted from the data, and the seven hyperfine components of the N2H + (1- 
0) spectra were fit simultaneously with Gaussians for emission and absorption (Caselli, Myers, 
& Thaddeus 1995). The properties of the N 2 H + (l-0) transition at the central position are 
shown in Table 1 and the spectrum at that position is shown in Figure 2. 

3. Dust Properties 

Thermal emission from dust can be used to determine such physical properties as tem- 
perature, density, and mass. However, the derivation of these quantities requires knowledge 
of the somewhat poorly known emission characteristics of the dust grains. In this section we 
discuss our assumptions and methods to derive the properties of TMC-1C from our sub-mm 
continuum observations. 
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3.1. Assumed Parameters 

The dust is assumed to emit as a modified blackbody, with emissivity parameter Q = 
(X/Xo)^ 13 . The flux coming from the dust at a particular wavelength is therefore given by 
(Mitchell et al. 2001) 

S x = nB x (T d )K xl 2m H N H (1) 

where 

2hc 2 1 

B\{T) = A 5 exp{hc/XkT) - 1 (2) 

and 

« A = K i300 (3) 



1300/iin, 

In Equation 1, 5^ is the flux per beam; Q is the solid angle of the beam; B\(Td) is the black 
body emission from the dust at temperature T<j; K1300 — 0.005 cm 2 g -1 is the emissivity of 
the dust grains at 1300 /im (Andre, Ward-Thompson, & Motte 1996; Preibisch et al. 1993); 
ijih is the mass of the hydrogen atom; \x = 2.33 is the mean molecular weight of interstellar 
material in a molecular cloud, and Nh is the column density of hydrogen nuclei. For a true 
blackbody (3 = 0; for amorphous, layerlattice material (3 ~ 1; for metals and crystalline 
dielectrics f3 ~ 2, cf (Henning, Michel, & Stognienko 1995). We assume the emission is 
optically thin at both 850 and 450 /im, which is justified for H2 column densities less than 
~ 10 24 (~ 500 magnitudes of visual extinction) (Zucconi, Walmsley, & Galli 2001). 



3.2. Derivation of Temperature 

Because the thermal emission from cold (~10 K) dust at 450 and 850 /im is both 
optically thin and near the peak of the modified blackbody spectrum, SCUBA is well-suited 
to observations of starless cores. A dust color temperature can be determined from the ratio 
of fluxes at two wavelengths. At temperatures near 10 K, the blackbody spectrum peaks 
around 300 fim, so for data taken at 450 and 850 /im the Rayleigh- Jeans approximation 
to the blackbody law cannot be used. The flux ratio therefore depends on the assumed 
emissivity spectral index and the temperature (Kramer et al. 2003) as: 

£450 = / 850\ 3+/3 exp(17K/r dust ) - 1 

5*850 

V4507 exp(32K/T, ust ) - 1 
where hc/Xk is 17 K and 32 K at 850 and 450 /tm respectively. 
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In order to derive the dust temperature in TMC-1C, we smoothed both SCUBA maps 
to the same 14" resolution. Because the 450 and 850 /iin SCUBA maps both have such high 
signal to noise, we were able to determine the temperature independently at each 14"xl4" 
pixel and make a detailed temperature map, assuming a constant value for (3 throughout the 
core 1 . By estimating a temperature at each point in TMC-1C from our observations, we can 
calculate self-consistent profiles for column density and mass. A map of the derived dust 
temperature is shown in Figure 3. The derived dust temperatures range from 6 K to 15 K. 



3.3. Calculation of Extinction 

Given the dust color temperature at each position, the column density can be calculated 
from the observed flux, and from this, the equivalent visual extinction can be estimated. 
Equation 1 can be rearranged to determine the column density from the observed flux and 
derived temperature. 

N„ = — (5) 

The equivalent extinction is given by 

E(B-V) A v 
Av = N « N H E(B-V) (6) 

A V = 5.3 x 10" 22 A^ (7) 

where Nh/ 'E(B-V) = 5.8 x 10 21 cm~ 2 is the conversion between column density of hydrogen 
nuclei (for our assumed gas to dust ratio) and the selective absorption, and Ry = A V /E(B — 
V) = 3.1 is the ratio of total-to-selective extinction (Mathis 1990; Bohlin, Savage, & Drake 
1978). Note that this value of Ry assumes a particular "color" for the absorbing dust, which 
may be different from core to core, or even within a single core. The extinction map created 
using these assumptions, and the value of k in Section 3.1, is shown in Figure 4. The derived 
extinction reaches a value of ~ 50 mag Ay at the center of the map. A map of the derived 
column density is shown in Figure 4. 



3.4. Derivation of the Emissivity Spectral Index 



Instead of assuming a constant value of (3 to calculate the temperature map as in Section 
3.2, one can alternatively calculate the value of the emissivity spectral index at each point in 



This restriction is relaxed in Section 3.4. 
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the map by choosing a constant value for the temperature. Using multi-transition spectral 
line data, Tafalla et al. (2002) find that the core L1544 has a constant gas temperature of 
approximately 10 K (with spatial resolution of ~ 0.03 pc), so experimenting here with cal- 
culations involving a constant dust temperature is justified. Because the emissivity spectral 
index appears in the equations used to derive the temperature and mass of a core, knowl- 
edge of its variability within a core is clearly valuable to the understanding of the physical 
properties of TMC-1C and cores in general. (3 depends on the dust grain size distribution, 
the composition of the mantle, and the surface area to volume ratio of the dust particles, all 
of which can vary throughout a core (Ossenkopf & Henning 1994). 

The derived emissivity spectral index map for T = 10 K at all radii is shown in Figure 
5. In the central 2 arcminutes (~0.08 pc) of the TMC-1C core, the value of (3 takes values 
in the physically plausible range 0.5 < (3 < 2.0. To self-consistently determine the column 
density, temperature, and emissivity spectral index of the dust simultaneously, we would 
need at least one more continuum map at another wavelength. For the remainder of this 
paper we will assume that (3 has a constant value and that the dust temperature in TMC-1C 
is not fixed, unless stated otherwise, although we understand that both quantities are likely 
to vary throughout the core. 



3.5. Calculation of Mass 



For optically thin emission, Equation 1 can be rearranged to convert a measurement of 
dust emission to mass, so that 



M = nfim H N H d 2 



S,d 2 



K\B\(Td) 



(8) 



where M is the mass of an emitting volume and d is the distance to it. Using the conversions 
given above and assuming a distance to Taurus of 140 pc, this equation is equivalent to 



M = 8.29 x 10~ 3 S 450 



exp 



32 K 



or 



M = 0.145 S S50 



exp 



17 K 



K 450 



0.026 cm 2 g^ 1 



K 850 



-1 



d 



140 pc 



0.01 cm 2 g- 1 



d 



140 pc 



(9) 



(10) 



The mass contained in each 14"xl4" pixel is determined by Equations 9 and 10. Figure 
6 shows the mass profile calculated by adding the mass in each pixel to annuli concentric 
around the peak column density. The mass in each pixel is determined using the 850 /iin 
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flux in each pixel and either the temperature assigned to each pixel shown in Figure 3 (the 
curve labeled "Actual Temperature Profile"), or a constant temperature assigned to every 
pixel. The total mass contained within a given radius of the TMC-1C column density peak 
can be converted into a density profile, assuming a three dimensional geometry, as discussed 
below. 



3.6. Temperature and Density Profiles 

Under the assumption of spherical symmetry, we can create three dimensional temper- 
ature and density profiles from the 850 and 450 /iin maps. To make the profiles, we break 
up the inner 2.2 arcminutes of TMC-1C into 14" wide concentric annuli, centered on the po- 
sition of the extinction peak. The outermost annulus only has flux contributed to it by the 
outermost spherical shell, and its temperature and density can be calculated as in Sections 
3.2 and 3.5. The next shell in recieves flux from the two outermost spherical shells. After 
subtracting the flux contributed from the outermost shell, whose temperature and mass has 
already been calculated, the remaining flux at 450 and 850 /im comes from the second spher- 
ical shell, allowing its temperature and density to be derived. In this manner we are able to 
work our way from the outside of TMC-1C in, finding the temperature and density of each 
spherical shell. The procedure is very similar to that described in David et al. (2001) for the 
Hydra A cluster of galaxies. 

In Section 3.5, the mass was calculated for each element in the map independently of 
its neighbors by assuming that a single temperature could be assigned to the entire column 
of material. In this section, instead of dividing the central 2.2 arcminutes of TMC-1C into 
isothermal cylinders, we break it up into isothermal spherical shells. The temperature and 
density profiles are fit by broken power laws. Each broken power law is described by five 
parameters: two coefficients (for normalization), two exponents (the slopes), and the break 
radius (where the two lines meet). The temperature and density profiles created in this way 
are shown in Figure 7. 

The derived temperature profile agrees well with the theoretical temperature profile of 
an externally heated Bonnor Ebert sphere presented in Gongalves, Galli, & Walmsley (2004), 
as well as the observed dust temperature profile of L183 (Pagani et al. 2004). Outside the 
break radius, the density profile of TMC-1C derived from the SCUBA data is similar to 
density profiles derived from N2H + (l-0) data and continuum dust emission of other starless 
cores (Table 2). Inside the break radius, TMC-1C has a considerably shallower density profile 
than these other starless cores. The starless cores to which we are comparing TMC-1C did 
not have their density profiles determined with the same concentric ring algorithm as used 



- 8- 



here, because researchers lacked the data to make the detailed temperature map that we 
were able to get from SCUBA. The power law slopes of the density profiles of these other 
cores can be thought of as isothermal approximations to the "true" density profile. For 
comparison, the powerlaw fit to the density profile of TMC-1C under the assumption of a 
constant T^ust = 10 K and using only the 850 /im flux is also shown in Figure 8, with the 
broken powerlaw fit shown in Table 2. 



The balance between outward turbulent and thermal pressures and inward gravitational 
pull determines the stability of a dense core. In this section we derive the virial mass from 
N2H + (l-0) data and show that for a wide range of density profiles and shapes the virial mass 
is significantly lower than the mass derived from dust emission. The N2H + (l-0) emission 
peaks in the same place that the dust column density peaks, so using this line to estimate 
the virial mass is justified. The N 2 H + (l-0) self-absorption profile shows evidence of infall 
motions, further suggesting the gravitational instability of TMC-1C (Schnee et al. 2005). 



The virial mass of the inner core of TMC-1C can be estimated from the observed line 
width of the isolated component of the N 2 H + (l-0) spectra. The total velocity dispersion has 
a non-thermal and a thermal component, given by: 



4. 



Virial Mass and Infall 



4.1. 



Calculation of the Virial Mass 



°TOT — NT ' T 



(11) 



The thermal velocity dispersion is given by the gas temperature by: 




kT 



(12) 



The line width depends on the velocity dispersion as: 



AV = V8 ln2a TO T 



(13) 



Where otot is the total velocity dispersion, <j NT is the non-thermal component of 
the velocity dispersion, o~t is the thermal component of the velocity dispersion, k is the 
Boltzmann constant, T is the gas temperature (assumed to be 10 K, (Tafalla et al. 2002)), 
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H-tracer is the mass of the molecule (N 2 H + or H 2 ), and AV is the observed line width. Note 
that the dust temperature is below 10 K throughout most of the TMC-1C core, so the virial 
mass calculated for a 10 K gas is an upper limit, if the gas and dust temperatures are coupled. 

The N 2 H + (l-0) spectra were annularly averaged with bins of 20", and the line width 
was then determined as a function of distance from the core's center. The N 2 H + (l-0) line 
width is ~ 0.25 km s" 1 , and is very nearly constant at all radii (as expected for a coherent 
core). The virial mass of an ellipsoid with a power law density profile determined in this 
way is 

= 5R (AvUH 2 )Y j_ 
vir G 8 In 2 a,a 2 [ } 

where Av to t is the quadratic sum of the 10 K thermal line width of H 2 and the non-thermal 
line width of N 2 H + (l-0), a± is a parameter than depends on the density profile, and a 2 is a 
parameter that depends on the ellipticity of the core. The virial mass is plotted as a function 
of radius in Figure 6, with the assumptions that the density profile is uniform (ai = 1) and 
the core is spherical (a 2 = 1). More sophisticated geometric assumptions are discussed in 
Sections 4.1.1 and 4.1.2. 



4-1.1. Dependence on the Density Profile 

We have shown in Section 3.6 that the density in TMC-1C is not uniform, and our 
estimate of the virial mass should take this into account. For a given line width, a more 
centrally condensed cloud will have a lower virial mass than a uniform density cloud (Bertoldi 
& McKee 1992). This can be parameterized for a power-law density distribution p(r) oc r~ k 

a - {1 ~ k/3) (15) 
ai ~ (1 - 2*75) (15) 

(Bertoldi & McKee 1992). 

The dependence of the virial mass on the density power law is shown in Figure 9. For 
r _1 and r~ 2 density distributions, the corrected virial masses are 90% and 60% of the value 
calculated by Equation 15 for a uniform density profile. Thus, the steeper the density profile, 
the greater the imbalance between the virial mass and the dust derived mass for TMC-1C, 
making it more unstable. 
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4-1.2. Dependence on the Shape of the Core 

Our estimate of the virial mass also depends on the shape of TMC-1C. It has been 
shown that many clouds are elliptical in projection, and usually prolate in three dimensions 
(Myers et al. 1991). Oblate clouds have slightly lower virial masses than a sphere whose 
radius is determined by the plane-of-sky dimension and a prolate cloud will have a larger 
virial mass than such a sphere. If the axial ratio, y, is less than 1 (an oblate cloud) then 

arcsin(l - y 2 ) 1 ' 2 
° 2 = (l- y 2)i/2 ( 16 ) 

If the axial ratio, y, is greater than 1 (a prolate cloud) then 

_ arcsinh(y 2 - 1) 1/2 

(Bertoldi & McKee 1992). The dependence of the virial mass on the axial ratio is shown 
in Figure 9. For axial ratios less than 10:1, the correction for the virial mass is less than a 
factor of ~ 3. The virial mass as a function of shape and density profile for TMC-1C arc 
shown in Figure 10, as are the masses derived from the dust emission for three representative 
values of (3. Notice that to make the virial mass equal to the dust-emission implied mass at 
0.06 pc in TMC-1C (see Figure 6), a prolate core would need to have an axial ratio of 20:1. 
This is ruled out by our data, except for the very unlikely case where we are viewing a cigar 
end-on. 



4.2. N 2 H+(l-0) Self-Absorption 

N2H + (l-0) is typically used as an optically-thin, high-density, low-depletion tracer of the 
kinematics of the interior of cores (Williams & Myers 1999; Caselli et al. 2002a,b). However, 
there is some evidence of N2H + (l-0) self-absorption in dense cores (Caselli et al. 2002b; 
Williams & Myers 1999). In another paper we will show that N 2 H + (l-0) is a better match 
to the dust emission in TMC-1C than most carbon bearing molecules, which tend to deplete 
onto dust (Schnee et al. 2005). In this paper we show evidence of redshifted self-absorption 
in N 2 H + (l-0) at several positions near the peak of the dust emission. 

To determine if the splitting of the N 2 H + (l-0) hyperfine components is due to self- 
absorption or two superimposed velocity components from different emission components 
along the same line of sight, the seven hyperfine components can be compared with each 
other to see if the location of the dip is constant amongst them. If this is the case, then 
kinematics are likely to be creating the non-gaussian spectra, if not (with the shift increasing 
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with the statistical weight of the component), then self-absorption is probably the cause. The 
dip in the TMC-1C N 2 H + (l-0) spectra does shift with statistical weight in a way consistent 
with self-absorption, so we believe that the double peaked spectra are not the result of two 
kinematically distinct clouds. While N 2 H + (l-0) self-absorption is uncommon enough to be 
interesting on its own, we use the N2H + (l-0) data in this paper to determine infall velocity 
of material close to the center of the TMC-1C core. 

The isolated component of each N 2 H + (l-0) spectrum (see Figure 2) was fit with two 
gaussians (one positive and one negative) to find the approximate velocities for the emission 
and absorption features, respectively. If the interior of the core is collapsing or flowing toward 
the center faster than the exterior, then radiative transfer models show that the spectra 
should have two peaks, with the blue peak higher than the red (Zhou, et al. 1993; Myers et 
al. 1996). Therefore, in an infalling region, the gaussian fit to emission will be bluer than the 
gaussian fit to absorption. In Figure 11 we show the histogram of the difference between the 
emission and absorption velocities for N2H + (l-0), where the absorption velocities are nearly 
all redshifted relative to the emission velocities, indicative of infall. Using the simple two 
layer radiative transfer model of a contracting cloud developed by Myers et al. (1996), we 
find that the infall velocity averaged over the central 0.05 pc of TMC-1C is 0.06 km s -1 . 

5. Discussion 

5.1. 450 and 850 /xm Emission Maps 

The 450 and 850 /im maps of TMC-1C look qualitatively different in the lower half of 
the maps. The 850 /im map has a peak at (-l',l'), though the 450 /im map has no peak 
there. The 450 /im map peaks at (0',-1.5'), where the 850 /im flux is comparatively weak. 
This suggests that the physical conditions in TMC-1C are different from the pre-protostellar 
cores L1544 and L1698B, in which the 450 and 850 /im fluxes are well correlated (Shirley et 
al. 2000). 

5.2. Temperature Distribution 

In order to compare the observed dust temperature in a starless core to theoretical 
models of externally heated cores, it is necessary to have high quality maps in at least two 
wavelengths. Our SCUBA maps of TMC-1C provide some of the best available submillimeter 
data of a starless core. In Figure 7 we show that the temperature profile is coldest at the 
center, with a temperature of ~ 6 K. The temperature rises by ~ 1 K out to 0.03 pc, after 
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which it rises more quickly, to a ~ 12 K at a radius of 0.08 pc. 

Our observed dust temperture profile is a close match to that shown in Figure 2 of 
Gongalves, Galli, & Walmsley (2004), in both temperature range and shape of the profile. 
Similarly, our dust temperature profile also matches that in Figure 4b of Evans et al. (2001). 
Both of these models simulate the temperature profile of externally heated starless cores 
with Bonnor-Ebert density distributions, which is a reasonable approximation of TMC-1C. 

5.3. Dependence on j3 

The parameters describing TMC-1C that we have derived above all assume a constant 
emissivity spectral index throughout the cloud, although a range of values of (5 have been 
tested, between 1.0 and 2.0. We have found that constant values of (3 smaller than 1.0, the 
mass derived from the temperature and the 450 /im flux is not consistent with the mass 
derived from the temperature and the 850 /im flux. In general, for fits with larger /3, the 
derived mass of the cloud increases and the temperature of the cloud decreases. For a low 
value of (3 of 1.0, the virial mass for the inner 0.06 pc of a uniform density sphere is still 2.5 
times less than the mass derived from dust emission (see Figure 6). So, although the dust 
derived mass varies strongly with f3, the conclusion that the virial mass is smaller than the 
dust derived mass of TMC-1C is robust (See Figure 10). 

5.4. Dust Emissivity 

The emissivity of dust has been measured with a combination of FIR flux maps and NIR 
extinction maps in the clouds IC 5146 and B68 (Kramer et al. 2003; Bianchi et al. 2003). 
Both studies find 850 /im fluxes consistent with k 850 = 0.01 cm 2 g _1 , with uncertainties in 
the range of 30% - 60%. TMC-1C is a starless core similar to B68 and those found in IC 
5146, so we estimate that uncertainties in our 450 and 850 /im calibrations (see section 2.1) 
and emissivities result in a derived mass that is accurate to within a factor of 2 for a given 
emissivity spectral index and dust to gas conversion. The masses we derive from the dust 
emission are greater than the virial mass by more than a factor of two, so the uncertainty 
in the dust emissivity can not by itself change the conclusion that TMC-1C should be in a 
state of collapse. However, a perverse conspiratorial combination of f3, axial ratio, and k v 
could still result in a dust mass lower than the N 2 H + (l-0) derived virial mass. Figure 10 
shows contours of the virial mass as a function of ellipticity and density profile, and also 
plots the mass derived from the 450 and 850 /iin fluxes for representative values of (3. It 
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is shown that for any density profile and any value of (3 tested, TMC-1C is gravitationally 
unstable for a wide range of axial ratios. If TMC-1C has a nearly spherical geometry, then 
the dust emissivity we use would have to be low by an unlikely factor of ~ 3-5 to make the 
virial mass and the dust derived mass approximately equal. 

5.5. Emissivity Spectral Index 

The emissivity spectral index, (3, can vary through a core, though it is often treated as a 
constant in papers on starless cores. By assuming a constant temperature (10 K) and using 
the 450 jum to 850 fim flux ratio, we have created a map of (3 (Figure 5), in the same way 
in which the temperature map 3 was created. The values of (3 are closer to zero (i.e. the 
dust emits more like a blackbody) towards the center of the core than towards the edges. 
This would be expected if the dust grains are larger at smaller radii in TMC-1C, which 
might be taken as evidence of grain growth (Testi, Natta, Shepherd, & Wilner 2003). On 
the other hand, it is also plausible that the emissivity spectral index does not change much 
throughout the core, and instead the temperature of the dust is lower in the center of the 
TMC-1C than towards the edges (Zucconi, Walmsley, & Galli 2001). This would be expected 
if interstellar radiation is able to heat the exterior of the core, but the extinction is too high 
in the center to allow for efficient heating there. Note the similarity in the dust temperature 
map (Figure 3) and the emissivity spectral index map (Figure 5). The emissivity spectral 
index for constant T = 10 K has values in the physically plausible range of 0.5 to 2.0 (center 
to edge, respectively), making it hard to call either the constant (3 or constant T assumption 
unreasonable. In order to solve for the dust column density, temperature, and emissivity 
spectral index we would need a third wavelength in our SED. 

5.6. Mass Estimates With and Without Maps at Multiple Wavelengths 

In most SCUBA observations, non-ideal weather degrades the quality of 450 /xm data 
significantly more than the 850 /im data. In such cases, the mass of a core is usually derived 
from just 850 /zm data by assuming that the core is isothermal and assuming a temperature 
and emissivity spectral index. In TMC-1C, because we can derive a self-consistent temper- 
ature map (assuming a value for f3) from 850 and 450 jum data, we can also derive the mass 
as a function of radius without assuming a constant temperature. 

Here we compare the mass that we derive from our 450 and 850 /im emission maps 
with what would ordinarily have been derived from a single wavelength map. The mass that 
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would have been derived from just our 850 /xm data and a constant temperature is shown in 
Figure 6 for temperatures of 5, 7, 10, 15, and 20 K. Note that the "true" cumulative mass 
profile lies somewhere between the 5 K and 10 K constant temperature mass profiles. To 
constrain this range further, we use the mass profile, M(r), derived from both the 450 and 
850 /iin fluxes along with the observed 850 /im cumulative flux profile, F(r), to determine 
a self-consistent temperature profile, T(r), which is a constant temperature within R. This 
profile answers the question "Given the total 850 /im flux contained within a cylinder of 
radius r and the mass contained within that radius, what isothermal temperature would 
one need to assume to make the mass and flux consistent with each other?" The result is 
shown in Figure 12. The values range from 7 K to match the inner 0.02 pc of a (3 = 2.0 
cloud to 11 K to match the mass out to 0.09 pc of a (3 = 1.0 cloud. This tight range 
on the allowed values for a uniform temperature dust model to match the mass derived 
from multiwavelength data demonstrates the difficulty of correctly guessing an appropriate 
temperature for an isothermal core and shows the importance of having multiple wavelength 
bands of observations to derive it directly. 

The range of isothermal temperatures typically assumed for cores which do not have mul- 
tiwavlength data available lies between 30 Kelvin for cores in Orion (Johnstone et al. 2001; 
Mitchell et al. 2001) and Ophiuchus (Johnstone et al. 2000) to 10 for cold pre-protostellar 
cores (Shirley et al. 2000). Our analysis of TMC-1C shows that the interior temperature of 
molecular cloud cores can be significantly colder than these assumed temperatures, as has 
been recently shown to be the case for other cores as well (Evans et al. 2001). The mass of 
gas and dust in TMC-1C calculated from the 850 /iin flux map for various isothermal cores 
is plotted in Figure 6. If we had assumed that the dust temperature was uniform at 20 K, 
then the mass derived from the dust emission would be smaller than the virial mass for any 
of the values of (3 that we considered. A 15 K cloud would have a dust mass nearly equal 
to the virial mass for 1.0 < (3 < 2.0, and a 10 K cloud would have a dust mass considerably 
higher than the virial mass for all values of (3. Because the inferred gravitational stability of 
a molecular cloud core is highly dependent on the assumed temperature, observations should 
be made at two or more wavelengths to confidently describe the state of a core. 

5.7. Observed and Model Profiles 

Though we have nowhere required that the resulting temperature, density, or flux profiles 
of TMC-1C be well fit by powerlaws, for both the temperature and density profiles, broken 
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power laws are good fits to the data 2 . In section 4.2 we found the infall velocity implied 
by the N 2 H + (l-0) spectra, without an attempt to explain the physical significance of the fit 
parameters. In this section we compare the predictions of several popular analytical star 
formation models with the data in order to predict the evolution of the TMC-1C core. Each 
model considered is fit to the density shown in Figure 7. The fit is done by minimizing the 
sum of the squared errors between the data values at each radius, and the predicted density 
at that radius. We have shown that in TMC-1C the inner density profile n(r) oc r~ as , the 
outer density profile is n(r) oc r -18 , and the infall speed is ~ 0.06 km s -1 . 

The inside out collapse model of a singular isothermal sphere (Shu 1977) predicts a 
broken power law density distribution with an r -2 profile outside the infall radius, and a 
profile inside the infall radius asymptotically approaching r~ L5 . Our observations also show 
a two-component density distribution, but the inner slope of the TMC-1C density profile is 
definitely shallower than predicted by the Shu model. To within the errors (see Table 2), 
the straight-forward inside-out collapse model can be ruled out by the density distribution. 
In addition, the infall velocity of the inside out collapse is predicted to be around the sound 
speed (a = 0.22 km s -1 for 10 K gas), while we have observed an infall speed of only 0.06 
km s _1 . 

The logotrope describes a pressure-truncated, self-gravitating sphere with the equation 
of state P/P c = 1 + A\n(p/p c ) (McLaughlin & Pudritz 1996). This model predicts a density 
profile that goes as r^ 1 in the outer parts of cores and much shallower near the center 
while the core is in equilibrium. Once the cloud begins to contract, an r -1,5 profile develops 
inside the expansion wave. The density profile in TMC-1C outer portions of TMC-1C is 
considerably steeper than that predicted by the logotropic model, and it is therefore also 
inconsistent. 

In many ambipolar diffusion models (Ciolek & Mouschovias 1994; Basu & Mouschovias 
1994), the density of an axisymmetric magnetically regulated cloud can be described by a 
broken power law. In these models, there is an inner region with a uniform density profile, 
surrounded by an envelope that has a power law slope ranging from r~ 2 immediately outside 
the uniform density central region and flattening to a mean ~ r~ 14 profile. This qualitatively 
agrees with our finding for TMC-1C, in that the inner region of the core has a flatter 
density profile than the outer region, though our derived inner powerlaw slope is significantly 
shallower than the ambipolar diffusion models predict. Certain ambipolar diffusion models 



2 The radius of the innermost point is taken to be one half of the HWHM of the 850 /im beam. Though the 
location of this point is somewhat arbitrary, its influence on the goodness of fit to the various star formation 
models is minimal. 
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developed to match the conditions (mass, density, rotation, magnetic braking) of the starless 
core L1544 do a slightly better job of describing TMC-1C than the untailored models. These 
models predict a steeper profile (~r~ L8 ) outside the break radius and a flattened interior 
profile that match a large portion of the TMC-1C density distribution, but fail to match the 
inner density point by wide margin. (Crutcher et al. 1994; Safier, McKee, & Stahler 1997). 
In the Crutcher model the infall speed is as low as 0.033 km s _1 at the boundary of the 
supercritical core, and as high as 0.133 km s -1 near the center (7x 10~ 4 pc). The infall speed 
in the central arcminute of TMC-1C, as calculated by the N 2 H + (l-0) spectra, is roughly 0.06 
km s _1 , in agreement with this model. 

The L arson- Penst on model describes the uniform collapse of an isothermal cloud. Like 
the expansion wave solution of Shu, the density in the outer regions of the cloud goes as 
r~ 2 , but it has a much flatter profile towards the center. This qualitatively agrees with the 
observed TMC-1C density profile. The infall velocity in the Larson-Penston model should 
be ~3.3a, where a is the sound speed (a = 0.22 km s _1 for 10 K gas). The TMC-1C infall 
velocities are far smaller than the ~ 0.7 km s^ 1 infall velocity predicted by the Larson- 
Penston model, and so it too is ruled out. 

The Bonnor-Ebert model describes a pressure bounded isothermal sphere, with a density 
profile that is non-singular at the origin (Ebert 1955; Bonnor 1956). Bonnor-Ebert spheres 
have a density profile that is close to r~ 2 at large radii, and flattens at smaller radii, similar 
to the Larson-Penston model (Harvey et al. 2001) and qualitatively in agreement with the 
TMC-1C density profile. The collapse of a Bonnor-Ebert sphere has been studied numerically 
by Foster & Chevalier (1993). Once the cloud has begun to collapse, their model predicts an 
expanding region of supersonic inflow, which asyptotically converges to the Larson-Penston 
flow. TMC-1C has no evidence of supersonic inflow, so this model fails to describe the core, 
though a Bonnor-Ebert sphere that is just beginning to collapse would not necessarily be 
expected to have supersonic infall motions. Furthermore, the density profile of the Bonnor- 
Ebert sphere is the best fit to the dust-derived density profile for TMC-1C. 

A summary of which models can be ruled out, or not, by their predicted density profiles 
and infall velocities is presented in Table 3 and Figure 13. 

6. Summary 

6.1. 450 and 850 /im data 

In sub-mm observations, 450 /im data are often degraded by weather far more than 850 
/im emission. In this paper we have shown that obtaining high quality 450 /im data, giving 
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us two wavelengths near the thermal peak of TMC-lC's modified blackbody spectrum, is 
critical to determining the temperature and mass distribution within a pre-stellar core. 

6.2. Low Temperature 

The temperature map shows that TMC-1C has a strong temperature gradient if (3 is 
roughly constant, with the inner regions of the core being colder than the outer by roughly 
a factor of two. This suggests that no point source has formed in TMC-1C to heat the cloud 
from the inside. The interior temperature of TMC-1C is significantly colder than 10 K for 
most values of f3, and rises to ~15 degrees at a radius of 0.1 pc. This temperature range is 
lower than typically assumed for star forming cores, and leads to a much higher dust and 
implied gas mass than would have been derived making a "typical" assumption of constant 
Tdust of 15 or 20 K. The temperatures derived for TMC-1C are consistent with theoretical 
models of pre-protostellar cores heated by an external attenuated radiation field (with or 
without cosmic rays), which predict cores with a ~7 K center heating up to ~14 K at the 
edges (Zucconi, Walmsley, & Galli 2001; Evans et al. 2001; Galli, Walmsley, & Gongalves 
2002). 

6.3. High Mass 

The virial mass of TMC-1C has been calculated from the line width of the N 2 H + (l-0) 
transition averaged over several annuli. We show that the total core mass implied by sub- 
mm dust emission is significantly higher than the virial mass for a wide range of density 
distributions and axial ratios. The virial mass at a radius of 0.06 pc is shown in Figure 10 
as a function of the exponent in the density power law and ratio of the axes. 

If instead of deriving the temperature from multiwavelength SCUBA data we had used 
an assumed temperature of 15 K, this conclusion would be far less certain. An assumed 
constant temperature of 20 K would have led us to conclude that the virial mass is larger 
than the mass determined from dust emission. 



6.4. The Effect of Varying Beta 

For values of the emissivity spectral index in the range from 1 to 2, the mass of the 
central 0.06 pc of TMC-1C varies from 9 M to 15 M , as determined by the same method 
used to create Figure 6. Choosing any value of (3 in this range does not change the conclusion 
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that the dust derived mass is much larger than the virial mass. Furthermore, changing the 
value of P does not significantly change the shape of the broken power laws that we fit 
to the temperature and density profile, and therefore does not change our fits to various 
star formation models. Using different values of P does change the absolute values of the 
temperature and density, however, as shown in Figure 7. If P is allowed to vary with radius 
and if the temperature of TMC-1C is held constant at 10 K throughout the inner 0.06 pc of 
the core, then the value of P needs to vary by over a factor of ~ 3 to account for observations, 
increasing towards the edges (see Figure 5). In reality, it is likely that P and T dust and N dusi 
all vary with radius, but with only two wavelengths observed, we are forced to hold one 
quantity fixed to calculate the other two. With observations at a third wavelength, this 
restriction can be relaxed, and all three quantities (P, Tdust, and Nd ust ) can be calculated 
simultaneously. 



6.5. Infall Motions Detected 

Self-absorbed N 2 H + (l-0) spectra show evidence of infall motions of roughly 0.06 km s' 1 
over a radius of ~0.05 pc. We believe that N 2 H + (l-0) is a good tracer of infall in TMC-1C 
because it is a high density tracer, its shape matches the dust emission near the column 
density peak, and because it shows no evidence of being significantly depleted. We therefore 
conclude that the interior of TMC-1C is flowing inwards at about 0.06 km s _1 . 



6.6. Star Formation Model 

TMC-1C is a starless core that is potentially beginning to collapse. Its virial mass is 
significantly lower than its dust derived mass, and N 2 H + (l-0) spectra show signs of sub- 
sonic infall. The derived density profile behaves qualitatively like that predicted by several 
star formation models in that it is consistent with an r" 2 power law outside of its break 
radius and is shallower interior to the break radius. However, none of the models tested here 
convincingly match both the density profile and infall velocity that we have measured (see 
Table 3 and Figure 13). 

A Bonnor-Ebert sphere provides the best fit to our dust-derived density profile, but 
even for that model the inner point is only barely within the error bars. Do our results 
mean that all extant detailed models of the collapse of a core into a star are "wrong" ? Not 
completely. The mismatch between theory and observations here is likely caused both by 
none of the models, on its own, being exactly right, and remaining assumptions necessary 
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in the data interpretation being imperfect. Perhaps incorporating elements of one theory 
into another (e.g. ambipolar diffusion in a Bonnor-Ebert sphere) will produce more realistic 
theories. And, perhaps adding (at least) another wavelength dust map will modify the (3, 
N^ ust , and T^ust distributions (mildly), or other spectral lines tracing material even closer to 
the dynamical center of cores will reveal (slightly) higher infall speeds. The excellent weather 
in Hawaii the day we observed TMC-1C has given us some of the tightest constraints yet on 
the physics of isolated star formation in cores. Now we need to rise to the challenge offered 
by these fine data. 

We would like to thank Paola Caselli, Neal Evans, Eric Keto, Charlie Lada, Phil Myers, 
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Fig. 1. — The 450/im (top) and 850/zm (bottom) emission maps of TMC-1C. The SCUBA 
beam size is displayed in each map. The coordinates of the (0,0) position are RA=4:41:38.8 
DEC=+25:59:42 (J2000). 



-24- 




Fig. 2. — N 2 H + (l-0) spectrum from the (0,0) position. Note the seven hyperfine components. 
The isolated component is shown at ~ -2.8 km s _1 . The velocity scale is only relative to the 
LSR for the central component. 
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Fig. 3. — Dust temperature derived from the 450 and 850/im SCUBA maps with the choice 
of p= 1.5. 
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Fig. 4. — Column density derived from the 450 and 850 /zm SCUBA maps with the choice 
of (3 — 1.5. The column density of particles is given by A^ pai . ticles = iV H /2.33, where iV H is 
derived as in Section 3.3. 
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Fig. 5. — Emissivity Spectral Index (j3) of the dust as determined by the 450 and 850 /j,m 
SCUBA maps with the assumption that the dust temperature is constant at 10K. 
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Observed and Virial Masses 
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Fig. 6. — The masses plotted here are those contained with a given radius of the central 
sub-mm extinction peak. The "virial mass" derived from the N 2 H + (l-0) line width is plotted 
(assuming a uniform density profile and constant temperatuer of 10 K), as is the mass derived 
from the 850 /im SCUBA map coupled with the "actual temperature profile," and the masses 
that would have been derived assuming various constant temperatures and the 850 /im data 
alone. 
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Fig. 7. — These plots show the dependence of density (top) and temperature (bottom) on 
radius from the center of TMC-1C, which is assumed to be a sphere. The lines show the 
best fit broken power law, as described in Section 3.6. It is assumed here that j3 = 1.5 
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Fig. 8. — These plots show the dependence of density (top) on radius from the center of 
TMC-1C, which is assumed to be an isothermal sphere. The density profile is shown in the 
top panel for comparison with Figure 7. The lines show the best fit broken power law, as 
described in Section 3.6. It is assumed here that (3 = 1.5 
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Fig. 9. — (Top) The dependence of the virial mass on the power law exponent of the density. 
The virial mass is smaller for more centrally condensed cores. (Bottom) The dependence of 
the virial mass on the axial ratio of an elliptical core, x < 1 is an oblate core, and x > 1 is a 
prolate core. Note that raising the virial mass by a factor of 5 requires a prolate shape with 
an axial ratio of 20:1 (see Section 4.1.2) 
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Virial Mass of an Ellipsoid with a Power Law Density Profile 
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Fig. 10. — The greyscale shows the dependence of the derived virial mass of TMC-1C 
contained within the central 0.06 pc of the core on the ellipticity and density power law 
exponent. The virial mass is calculated with the N 2 H + (l-0) linewidth as in Section 4.1. The 
highlighted contours show the mass of TMC-1C contained within the central 0.06 pc of the 
core derived from dust emission with /3 taking the values 2.0 (dashed), 1.5 (dotted), and 1.0 
(solid). The solid horizontal line shows the approximate prolate axial ratio of TMC-1C as 
seen in projection. The solid vertical line shows the best fit single power law density profile 
exponent. 
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Fig. 11. — The value of the difference between the velocities of the emission and absorption 
components of the N 2 H + (l-0) spectra, fit by two gaussians (see Section 4.2). 



-34- 




Fig. 12. — This plot shows the temperature that one would need to assume is constant within 
radius R in order to correctly derive the mass enclosed within that radius from just the total 
850 micron flux within that radius for various values of the emissivity spectral index. 
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Fig. 13. — Best fit density profiles for TMC-1C, for various models. 
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Table I. N 2 H+(l-0) Spectrum at (0,0) a of TMC-1C 



Transition 


Frequency 


Signal / Noise 


Noise 


Vls.r 


Line Width 




GHz 




Kelvin 


km s _1 


km s _1 


N 2 H+(l-0) 


93.176 


9.8 


0.145 


-2.787 


0.246 



a (RA,DEC) (J2000) at (0,0) 04:41:38.8 +25:59:42 
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Table 2. Density Profiles of Starless Cores 



Observation Type 


Inner Exponent 


Outer Exponent 


Break Radius (pc) 


Reference 


N2H+(l-0) starless cores 


1.2 


2 


0.03 


1 


L1696 sub-millimeter 


1.3 


2 


0.02 


2 


L1689B sub-millimeter 


1.0 - 1.4 


2 


0.02 


3 


Sub-millimeter starless cores 


1.25 


2 


0.02 


4 


TMC-1C sub-millimeter (3D) 


0.8 ± 0.1 


1.8 ± 0.1 


0.02 




TMC-1C sub-millimeter (3D, constant T^ ust ) 


0.6 


1.1 


0.02 





References. — (1) (Caselli et al. 2002b); (2) (Ward-Thompson, Motte, & Andre 1999); (3) (Andre, Ward-Thompson, & 
Motte 1996); (4) (Ward-Thompson ct al. 1994) 
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Table 3. Star Formation Models 



Model 


Density Profile 


Infall Velocity 


Bonnor Ebcrt Sphere 


not ruled out 


not ruled out 


Inside Out Collapse 


ruled out 


ruled out 


Larson Penston Model 


ruled out 


ruled out 


Logotrope 


ruled out 


not ruled out 


Ambipolar Diffusion 


ruled out 


not ruled out 



